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"Shear  Dispersion"  in 
Time- Varying  Flows 


W.C.  Thacker 


Abstract.  Shear  dispersion  in  time-varying  flows  is  shown  to  be  parameteriz- 
able  as  effective  diffusion  with  a  time-varying  diffusivity.  This  parameteriza- 
tion is  appropriate  only  for  times  long  in  comparison  with  the  cross-shear 
mixing  time.  High  frequency  variations  of  the  flow  do  not  contribute  substan- 
tially to  the  mixing.  Low  frequency  variations  of  the  flow  can  be  accounted  for 
by  simply  allowing  the  parameters  in  the  expression  for  the  effective  diffusivity 
to  take  on  a  time  dependence. 


1.  Introduction 

Previous  discussions  (Bowden,  1965;  Schonfeld,  1961;  Okubo,  1967)  of 
harmonically  varying  shear  have  concluded  that  high  frequency  time  variations 
do  not  contribute  substantially  to  the  mixing.  Here  it  is  shown  that,  further- 
more, a  time-dependent  diffusivity  can  be  introduced  to  parameterize  time- 
dependent  shear  dispersion.  A  simple  solvable  model  (Thacker,  1976)  is  used  to 
obtain  closed  form  expressions  for  shear  dispersion  with  arbitrary  time  depend- 
ence. 

2.  Closed  Form  Expressions 

Shear  dispersion  is  mixing  due  to  the  combined  effects  of  velocity  shear  and 
cross-shear  mixing.  A  simple  model,  which  views  the  shear  as  a  flow  consisting 
of  two  layers  of  fluid  moving  in  opposite  directions  with  velocities  ±  //and  with 
mixing  of  contaminant  between  the  layers  at  a  rate  a,  can  be  expressed  by  the 
equations 
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where  C]  (x,t)  andC2  (x,t)  are  distributions  of  contaminant  in  the  two  layers,  and 
■Si  (x,t)  and  S2  (x,t)  are  sources  of  contaminant.  If//  and  a  do  not  vary  in  time, 
thenC  =  V3  (Ci  +  C2)  is  well  approximated  by  the  solution  of  a  diffusion  equation 
with  effective  diffusivity  given  by  K*=//2/2a  (Thacker,  1976).  It  will  be 
shown  that  the  same  is  true  even  if//  and  a  depend  upon  time. 

For  the  case  of  time-independent  flow,  where  //  and  a  are  not  functions  of 
time,  closed  form  solutions  for  C  =  1/2(C)  +  C2)  and  AC  =  1/2(C]  —  C2)  can  be 
obtained  in  terms  of  a  Green's  function, 
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0;  otherwise. 
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The  matrix  notation, 
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allows  the  solution  to  ( 1)  for  arbitrary  initial  conditions C(x,t  =  0)  and  sources  of 
contaminant  S(x,t)  to  be  expressed  compactly  as 


C{x,t)=\X  dx    {dt'%>a(x-x',t-t')S(x',t')+  |    * dx,c6UM  (x-x',t)C(x',t=0).      (4) 
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Equation  (4)  can  be  shown  to  be  valid  even  when  //  and  a  are  functions  of 
time.   In  this  case,  the  matrix  Green's  function  is  given  by 
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The  matrices  ^w.,«.  are  the  matrix  Green's  functions  that  govern  the  solution  in 
/th  time  interval',  from  (/'  —  l)t /Nto  it/N  in  which  the  parameters  have  the 
approximately  constant  values  a,  and//,.  It  is  important  that  the  matrix  Green's 


functions  ^„.,u.  are  arranged  for  multiplication  with  the  index  i  decreasing  to 
the  right  because  they  do  not  commute.  Thus,  ^uWMt)  1S  the  N-fold  convolu- 
tion of  the  Green's  function  for  each  of  N  time  intervals  composing  the  total 
interval  from  0  to  t.  In  the  limit,  N  goes  to  infinity,  so  that  this  is  an  infinite 
convolution. 

This  infinite  convolution  is  difficult  to  evaluate.  However,  it  is  possible  to 
obtain  closed  form  expressions  for  the  moments  (X")  of  the  mean  contaminant 
distribution C(x,t)  for  arbitrary  functions  cx(t)  and //(/).  The  Fourier  transform  of 
C(x,t)  is  the  generating  function  for  these  moments.  Equation  (4)  expresses 
C(.v./)  in  terms  of  convolution  of  initial  conditions  and  sources  with  ^M(o,a«)> 
which  in  turn  is  an  infinite  convolution.  Since  the  Fourier  transform  of  a 
convolution  is  a  product,  the  problem  is  reduced  to  an  easier  one  of  evaluating  an 
infinite  product. 


The  Fourier  transform  of  ^uu),  a(t)  is  given  by 
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where  Q,i=(Ui  k      a,-  )*.  The  moments  of  the  distribution  can  be  obtained  by 
multiplying  moment  matrices,  which  are  generated  according  to 
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k  =  0 


with  the  Fourier  transform  of  the  initial  distribution,  if  there  are  no  sources. 
The  first  few  moment  matrices  are  found  to  be 
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Consider  the  case  of  initially  localized  distributions,  C(x,t  =  0)  =  8(.v)  and 
AC(.v,/  =  0)  =  /38(a),  I  )8  I  <  1.  If  /3  =  0,  the  contaminant  is  initially  distri- 
buted equally  between  the  two  layers;  if /3=  ±  1,  it  lies  entirely  within  one  of  the 
two  layers.  Also,  suppose  that  there  are  no  continuous  sources  of  contaminant 
S(x,t)  —  AS(x,t)  =  0.  Then  the  first  few  moments  of  C(x,t), 
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are  given  by 
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These  moments  can  be  evaluated  for  arbitrary  time  dependence  of  the  shear 
and  mixing.  The  zeroth  moment  indicates  that  the  total  amount  of  contaminant 
is  constant  in  time,  regardless  of  variations  of  the  flow.  The  first  moment 
indicates  that  there  is  a  horizontal  asymmetry  only  if  initially  there  is  a  vertical 
asymmetry.  However,  the  magnitude  of  the  horizontal  asymmetry  depends 
upon  the  time  dependence  of  the  flow  parameters,  //  and  a.  The  second  moment 
measures  the  square  of  the  width  of  the  distribution. 


3.     Parameterization  as  Time-Dependent  Diffusion 


A  diffusion  process  with  time-dependent  diffusivity  K*(t)  would  yield  the 
moments 
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Comparison  of  the  second  moment  in  ( 10)  with  that  in  ( 1 1),  shows  that  in  some 
sense  the  expression 
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must  approach  K*  (t)  if  the  time-dependent  shear  effect  is  to  behave  like 
time-dependent  diffusion.  This  becomes  more  clear  after  consideration  of  the 
following  cases  for  different  time  dependencies  of//  and  a. 


Case  One.  In  the  case  of  no  time  dependence,  where//  and  a  are  constant, 
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The  time  dependence  of /  should  not  be  interpreted  as  time-dependent  diffu- 
sion. Shear  dispersion,  lor  this  case,  can  be  parameterized  as  diffusion  only  for 
times  long  in  comparison  with  the  cross-shear  mixing  time;  i.e.,  at  »1.  In  this 


limit  K*=f=—,  (x2)=2(u2/2a)t,  and  (x2)»(x)2  (Thacker,   1976). 


Case   Two.    In    the   case   of    harmonically   varying    shear,    with  u{t)~- 
U),  cos(2777 fT) ,  and  constant  a. 


"A 


(*)  = 


(*2)  = 


B  —1 1 

P  2a1 


cos 


2737 

r 


7T       .      277/  \ . 


7*'  i  + 

4a 


I 
a/ 


IL  sin  — '  +  —  (■ 
47rr      '   T        4a/ 1 


cos 


477/ 

T 


277/ 

e  ic"\  cos  — 

r 


77       .      277/ 

a.T         T 


R\R        > 


/(')    -  7T-  COS^  — 

2a  T 


TT 


277/ 

1  +  —  tan  — 
aT         T 


277/    _2a?| 

sec  ? 

r 


— V2"'  i? 


A' 


(aT) 


(aT)  +  77J 


(14) 


y 


Each  expression  contains  a  factor  R,  which  indicates  that  high  frequency  shear 
does  not  contribute  to  the  dispersion,  since R—>0  as  aT^O.  This  result  was  also 
found  by  Schbnfeld  (1961)  for  this  model  and  by  Okubo  (1967)  for  another 
harmonic  shear  flow.  The  expressions  that  should  be  compared  with  time- 
dependent  diffusion  are  those  valid  in  the  long  time  limit  at»l.  With  R=  1, 
these  expressions  are 
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Thus,  low  frequency  harmonic  shear  can  be  parameterized  as  time-dependent 
diffusion  with  K*  =  \u(t)fl2a. 

Schonfeld  (1961),  Bowden  ( 1965),  and  Okubo  (1967)  have  attempted  to 
describe  harmonic  shear  as  time-independent  diffusion.  They  obtained  a  diffu- 
sion coefficient  given  by  ( \JT)\TK*  (t)  dt,  the  time  average  of  A'*  (/)  over  one 
cycle  of  the  shear.  This  yields  the  same  result  for  be  \  in  the  diffusion  limit,  but 
does  not  allow  for  generalization  from  harmonic  shear  to  arbitrary  shear. 

Case  Three.  The  question  remains  of  how  the  diffusivity  should  be 
parameterized  for  the  general  case  in  which  the  shear  and  mixing  vary  arbitrarily 
in  time.   The  answer  should  be  obtainable  by  analyzing  equation  (12).   This 


analysis  is  difficult  because  at  any  point  in  time  there  are  three  time  scales  that 
must  be  considered:  the  inverse  of  the  mixing  rate,  /]  =a~1;  the  time  over  which 
the  mixing  varies,  t2=a/(da/dt);  and  the  time  over  which  the  shear  varies, 
t^—u/(du/dt).  On  the  basis  of  the  results  for  cases  one  and  two  and  the 
discussion  of  Thacker  ( 1976),  it  seems  reasonable  that  the  result  of  the  analysis 
for  the  qeneral  case  should  be 

*•-££  ,16, 

2a(t) 

for  t2»t\  and  /,;»/] . 

4.     Discussion 

The  result  presented  here  is  that  shear  dispersion  can  be  parameterized  as 
enhanced  diffusion  along  the  shear  even  for  time-dependent  flows,  so  long  as 
sufficient  time  has  elapsed  for  the  cross-shear  mixing  to  take  place.  This  is  an 
extension  of  a  previous  time-independent  analysis  (Thacker,  1976).  All  of  the 
conclusions  from  that  analysis  should  also  apply  for  the  time-dependent  case. 
First,  the  diffusion  parameterization  of  the  mixing  should  be  used  when  the 
details  of  the  shear  are  not  resolved.  This  is  equivalent  to  ignoring  variations  in 
times  shorter  than  the  time  required  to  mix  contaminant  across  the  shear  and  to 
ignoring  variations  along  the  shear  in  distances  shorter  than  the  distance  over 
which  the  shear  stretches  the  contaminant  in  the  cross-shear  mixing  time. 
Second,  the  effect  of  turbulent  mixing  can  be  accounted  for  simply  by  adding  the 
turbulent  eddy  diffusivity  to  the  diffusivity  that  parameterizes  the  shear  disper- 
sion. Finally,  a  more  sophisticated  description  of  the  shear  flow,  such  as  that 
discussed  by  Bowden  (1965),  would  lead  to  a  more  complicated,  and  perhaps 
more  accurate,  prescription  for  relating  the  effective  diffusivity  parameter  to  the 
shear  flow.  However,  the  result  found  here  should  also  hold  for  more  sophisti- 
cated descriptions  of  the  shear  flow.  The  time-dependent  diffusivity  can  be 
obtained  from  the  expression  for  the  time-independent  diffusivity  simply  by 
allowing  the  flow  parameters  in  that  expression  to  vary  in  time. 


5.      References 

Bowden,  K.  F. ,  1965.  Horizontal  mixing  in  the  sea  due  to  a  shearing  current. 

J.  Fluid  Me ch..  21:83-95. 
Okubo,  A.,   1967.   The  effect  of  shear  in  an  oscillatory  current  on  horizontal 
diffusion   from   an    instantaneous   source.    Int.  J.    Oceanol.    Limnol.,    1: 
194-204. 
Schonfeld,  J.  C,    1961.   The  mechanism  of  longitudinal  diffusion  in  a  tidal 

river.  Bull.  Int.  Assoc.  Sci.  Hydro/.    (Louvain,  Belgium),  6(1). 
Thacker,   W.   C,    1976.    A  solvable  model  of  "shear  dispersion."  J.   Phys. 
Oceano^r.,  6:66-75. 


ft  U.S.  GOVERNMENT  PRINTING  OFFICE:  1977— 777-067/1218  REGION  NO.  8 


LABOR  AT  DRIES 


The  mission  of  the  Environmental  Research  Laboratories  (ERL)  is  to  conduct  an  integrated  program  of  fundamental 
research,  related  technology  development,  and  services  to  improve  understanding  and  prediction  of  the  geophysical 
environment  comprising  the  oceans  and  inland  waters,  the  lower  and  upper  atmosphere,  the  space  environment,  and  the 
Earth.  The  following  participate  in  the  ERL  missions: 


MESA         Marine  EcoSystems  Analysis  Program.  Plans, 
directs,  and  coordinates  the  regional  projects 
of  NOAA  and  other  federal  agencies  to 
assess  the  effect  of  ocean  dumping,  municipal 
and  industrial  waste  discharge,  deep  ocean 
mining,  and  similar  activities  on  marine 
ecosystems. 

OCSEA       Outer  Continental  Shell  Environmental 

Assessment  Program  Office.  Plans  and  directs 
research  studies  supporting  the  assessment 
of  the  primary  environmental  impact  of  energy 
development  along  the  outer  continental  shelf 
of  Alaska;  coordinates  related  research  activities 
of  federal,  state,  and  private  institutions. 

WM  Weather  Modification  Program  Office.  Plans, 

directs,  and  coordinates  research  within  ERL 
relating  to  precipitation  enhancement  and 
mitigation  of  severe  storms.  Its  National 
Hurricane  and  Experimental  Meteorology 
Laboratory  (NHEML)  studies  hurricane  and 
tropical  cumulus  systems  to  experiment  with 
methods  for  their  beneficial  modification  and 
to  develop  techniques  for  better  forecasting 
of  tropical  weather.  The  Research  Facilities 
Center  (RFC)  maintains  and  operates 
aircraft  and  aircraft  instrumentation  for 
research  programs  of  ERL  and  other  govern- 
ment agencies. 

AOML         Atlantic  Oceanographic  and  Meteorological 
Laboratories.  Studies  the  physical,  chemical, 
and  geological  characteristics  and  processes 
of  the  ocean  waters,  the  sea  floor,  and  the 
atmosphere  above  the  ocean. 

PMEL  Pacific  Marine  Environmental  Laboratory. 

Monitors  and  predicts  the  physical  and 
biological  effects  of  mans  activities  on 
Pacific  Coast  estuarine,  coastal,  deep-ocean, 
and  near-shore  marine  environments. 

GLERL        Great  Lakes  Environmental  Research  Labora- 
tory. Studies  hydrology,  waves,  currents,  lake 
levels,  biological  and  chemical  processes, 
and  lake-air  interaction  in  the  Great  Lakes  and 
their  watersheds;  forecasts  lake  ice  conditions. 


GFDL      Geophysical  Fluid  Dynamics  Laboratory. 
Studies  the  dynamics  of  geophysical  fluid 
systems  (the  atmosphere,  the  hydrosphere, 
and  the  cryosphere)  through  theoretical 
analysis  and  numerical  simulation  using  power- 
ful, high-speed  digital  computers. 

APCL       Atmospheric  Physics  and  Chemistry  Labora- 
tory. Studies  cloud  and  precipitation  physics, 
chemical  and  particulate  composition  of  the 
atmosphere,  atmospheric  electricity,  and 
atmospheric  heat  transfer,  with  focus  on 
developing  methods  of  beneficial  weather 
modification. 

NSSL       National  Severe  Storms  Laboratory.  Studies 
severe-storm  circulation  and  dynamics,  and 
develops  techniques  to  detect  and  predict 
tornadoes,  thunderstorms,  and  squall  lines. 

WPL         Wave  Propagation  Laboratory.  Studies  the 
propagation  of  sound  waves  and  electro- 
magnetic waves  at  millimeter,  infrared,  and 
optical  frequencies  to  develop  new  methods 
for  remote  measuring  of  the  geophysical 
environment. 

ARL         Air  Resources  Laboratories.  Studies  the 

diffusion,  transport,  and  dissipation  of  atmos- 
pheric pollutants;  develops  methods  of 
predicting  and  controlling  atmospheric  pollu- 
tion; monitors  the  global  physical  environment 
to  detect  climatic  change. 

AL  Aeronomy  Laboratory.  Studies  the  physical 

and  chemical  processes  of  the  stratosphere, 
ionosphere,  and  exosphere  of  the  Earth  and 
other  planets,  and  their  effect  on  high-altitude 
meteorological  phenomena. 

SEL  Space  Environment  Laboratory.  Studies 

solar-terrestrial  physics  (interplanetary,  mag- 
netospheric,  and  ionospheric);  develops  tech- 
niques for  forecasting  solar  disturbances; 
provides  real-time  monitoring  and  forecasting 
of  the  space  environment. 


U.S.  DEPARTMENT  OF  COMMERCE 
National  Oceanic  and  Atmospheric  Administration 

BOULDER,  COLORADO  80302 


PENN  STATE  UNIVERSITY  LIBRARIES 


II 


AODDDTZOEIAST 


